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Abstract 

We discuss reactor measurement of #12 which has a potential of reaching the ultimate sensi- 
tivity which surpasses all the methods so far proposed. The key is to place a detector at an 
appropriate baseline distance from the reactor neutrino source to have an oscillation maximum at 
around a peak energy of the event spectrum in the absence of oscillation. By a detailed statistical 
analysis the optimal distance is estimated to be ~ (50 — 70) kmx[8 x 10 _5 eV 2 /Am| 1 ], which is 
determined by maximizing the oscillation effect in the event number distribution and minimizing 
geo-neutrino background contamination. To estimate possible uncertainty caused by surrounding 
nuclear reactors in distance of ~ 100 km, we examine a concrete example of a detector located at 
Mt. Komagatake, 54 km away from the Kashiwazaki-Kariwa nuclear power plant in Japan, the 
most powerful reactor complex in the world. The effect turns out to be small. Under a reasonable 
assumption of systematic error of 4% in the experiment, we find that sin 2 9\2 can be determined 
to the accuracy of ~ 2% (~ 3%), at 68.27% CL for 1 degree of freedom, for 60 GW t h-kton-yr (20 
GWth - kton-yr) operation. We also discuss implications of such an accurate measurement of 6\2- 
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I. INTRODUCTION 



Reactor neutrino experiments have been the source of rich physics informations The 
very existence of (anti-)neutrino was proven by the memorial reactor neutrino experiment 
by Reines and Cowan |2j. In recent years, reactor neutrino experiments have been play- 
ing greater role than in any other era of neutrino physics. Following SNO 0, 0] who has 
confirmed by the in situ measurement that the solar neutrino deficit is indeed caused by neu- 
trino flavor transformation, KamLAND 0, 0] observed clear deficit of reactor neutrinos and 
pinned down a unique parameter region, the large-mixing-angle (LMA) Mikheyev-Smirnov- 
Wolfenstein (MSW) solution 0. Thus, the solar neutrino problem which lasted nearly 40 
years has now been solved. At the same time, the basic structure of the (1-2) sector of the 
lepton flavor mixing matrix, the Maki-Nakagawa-Sakata (MNS) matrix was determined. 

Reactor neutrino experiments also play (and will play) crucial role to explore the (1-3) 
sector of the MNS matrix, the unique sector whose structure is not yet determined. The 
CHOOZ and the Palo Verde experiments were able to place stringent constraints on 9 13 
0, [lij • It was recognized that reactor measurement of #13 has very special characteristics 
as a pure measurement of #13, whose prop erty should play an important part in solving 



the problem of parameter degeneracy [111 ]. It was also stressed that the role played by 
reactor experiments is complementary to those by long-baseline (LBL) accelerator neutrino 
experiments. The proposal, sometime after an earlier Russian project [l2|, was followed by a 
spur of experimental projects over the globe which are now summarized in the White Paper 
Report Q. 

Are those described above all that can be done by reactor neutrino experiments? We 
answer the question in the negative by proposing dedicated reactor neutrino experiments 
for precise measurement of (9 12 in this paper. We will show that the accuracy attainable 
by such reactor experiments for 9 12 can reach to ~ 2% for sin 2 9 12 for 1 degree of freedom 
(DOF) under a reasonable assumption of systematic errors. It surpasses those expected 
by the other methods so far proposed, e.g., by combining KamLAND with accurate 7 Be 
and pp solar neutrino measurement [14J . Throughout this paper, we will demonstrate these 
statements by careful treatment of the optimization of the baseline distance as well as the 
effects of geophysical neutrinos (hereafter, geo-neutrinos) . 

The key idea for such an enormous sensitivity is to place the detector at the appropriate 
distance to see the effect of oscillation at the place where it is maximal. A zeroth order 
estimation of the distance ^| is given byLoM = 2nE^ Qak / Am^ ~ 63 km, where Am^ = 
7.9 x 10~ 5 eV 2 is the current best fit value 0, and E? cak = 4 MeV is a typical neutrino energy 
where event rate has a maximum in the absence of oscillation. But, due to the fact that the 
reactor neutrino energy spectrum is rather broad and the 1/L 2 dependence of the intensity 
of neutrino flux, the optimal distances turned out to be shorter and spread over a range of 
50 to 70 km, adding more variety for the site selection of the detector. We also discuss the 
geo-neutrino background contamination and its relevance in determining the best position 
to place the detector. 

It seems that now a coherent view of how to determine accurately the "large" mixing an- 
gles has emerged. Namely, disappearance measurement at around the oscillation maximum 
gives the highest sensitivities for both of the large angles, 9 12 and #23 . It is well recognized 
that the most accurate way of determining sin 2 2^23 will be achieved by the next generation 
LBL accelerator neutrino experiments by using their disappearance mode [la . 17 . E^ |. 



In particular, in the JPARC-SK experiment, the accuracy of determination of sin 2# 2 3 is 
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expected to reach down to ~ 1% at 90% CL [3, 3 • 

What is the scientific merit of such precise measurement of #12? We will discuss this 
question in length in Sec. IVIII1 but here we make only two remarks, one from particle 
physics and the other from solar physics point of view. Current understanding of nonzero 
neutrino mass involves, most probably, the existence of a large energy scale where leptons and 
quarks are unified j^, whose simplest model realization is the seesaw model Various 
models so far proposed try to relate quark and lepton mass matrix based on the underlying 
philosophy. For a review see e.g., 2^\. We will probably need detailed informations of lepton 



mixing matrix which are comparable with that of quarks when we want to test the theory 
for lepton and quark mixings, which should emerge someday in the future. If the solar angle 
9 12 is in fact complementary to the Cabibbo angle |23|, there is an immense request for 
accurate determination of 9\i ■ 

What is the solar physics implication of accurate measurement of 9yp- The solar neutrino 
flux measured on the earth inevitably contains the effect of neutrino flavor transformation, 
and therefore the precision of its determination is affected by the uncertainty of 9\i- We 
thus believe that accurate determination of 9\2 has a great merit for precise observational 
solar astrophysics in which an accurate measurement of the infant flux produced at the solar 
core is required. 

In Sec. |Tl]we discuss the basic principle and power of the method of tuning the baseline 
distance for accurate determination of 6*12 by reactor neutrinos. In Sec. IIHI we discuss 
the requirements and possible obstacle for such reactor # 12 experiments. They include the 
systematic errors, a possible uncertainty due to geo-neutrino background and the unknown 
value of 9i 3 . In Sec. HVJwe describe the analysis procedure and give a quantitative estimate 
of the optimal baseline distances by fully taking into account the geo-neutrino background. 
We use a specific setting of the detector to estimate effects of the other surrounding reactors. 
In Sec. |?lwe carry out the sensitivity analysis. In Sec. IVII we examine the stability of our 
results against changes of our statistical procedure. In Sec. IVHI we compare the sensitivity 
of reactor #12 experiments with the one that will be reached by KamLAND-solar combined 
method. In Sec. IVIIII we discuss possible physics implications of a precise measurement of 

#12- 



II. ACCURATE MEASUREMENT OF 6 12 BY REACTORS BY TUNING BASE- 
LINE DISTANCE 

We discuss in this section the method for accurate determination of #12 by choosing an 
appropriate distance to the detector from a principal (nearest) nuclear reactor. We start by 
giving a pedagogical self-contained description of the basic principle of the measurement. 
We then contrast the sensitivity achievable by our method with that of KamLAND j^, 0| 
to illuminate the power of the method we propose by explicitly showing improvement over 
the marvelous experiment using the same neutrino flux from the reactors located all over 
Japan. 

For ease in frequent reference to a dedicated liquid scintillator detector at the distances 
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of 50-70 km, we use in this paper the acronym "SADO" as an abbreviation, 1 

SADO = Several-tens of km Antineutrino DetectOr. (1) 

A. Method for precise measurement of 9\2 

The secret behind the potentially powerful method is to place a detector to the location 
where the oscillation effect has its maximum in the observable quantity, i.e., the number of 
events (not in the probability alone). The optimal distances for a given value of Am^ are 

-l 

km, (2) 

as we will find in Sec. II VI after establishing our statistical method (See Fig. El) The optimal 
distance coincides approximately, but not exactly, with the first oscillation maximum (the 
one encountered first by neutrinos as they travel) for the peak energy of reactor neutrino 
event distribution. Moreover, such a distance must be established by taking into account 
background contribution from geo-neutrinos, as we will see in Sec. IIHI 

In FigGJwe show the range of survival probabilities as a function of neutrino energy in 
a form of band which is spanned by P{p e v e ) obtained by varying sin 2 #12 between 0.25 
(upper end) and 0.35 (lower end). The probabilities are computed using the current best fit 
value of Am^ and with the two distances L = 50 (red) and 60 km (blue). 

We observe that the oscillation maximum (minimum in the survival probability) occurs at 
E v ~ (3.0-4.0) MeV around the peak energy where the event rate is maximal in the absence 
of oscillation. The most important feature for us is that the depth of the minimum in the 
survival probability is the most sensitive place to the variation of sin 2 #12 in the entire figure. 
Nothing but that property is the key to the enormous accuracy of reactor #12 measurement 
to be explored in this paper. 

We should note that there were some previous attempts along the similar line of thought. 
To our knowledge, dedicated reactor experiment for # 12 at the oscillation maximum has first 
been considered by Takasaki who did a back of the envelope estimation of the sensiti vity 
assuming an appropriate site at ~70 km away from the Fukushima nuclear reactors [25]. 
A statistical treatment for estimating the accuracy of 6* 12 determination was presented by 
Bandyopadhyay et al. Jl5| which resulted in about 10% error (99% CL, 2 DOF) in sin 2 #i 2 . 
It is about factor of 2 larger than ours; See Fig. 0] They also did not address the optimiza- 
tion problem of the sensitivity with respect to the baseline distance. A different strategy 
which utilizes not only the first oscillation maximum but also the subsequent minimum was 
proposed by Bouchiat |26| . 

B. Setting 

Throughout this paper we analyze simultaneously the two settings: 



"SADO" was originally created as an acronym for a detector in Sado island in Niigata, Japan |24] , which 
is located at about 71 km from the Kashiwazaki-Kariwa nuclear reactor complex, whose distance roughly 
corresponds to the oscillation maximum for the old best fit value of Ato^ = 7.2 x 10~ 5 eV 2 0. 
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FIG. 1: Ranges of the survival probabilities P{p e — > P e ) of antineutrinos from the nuclear reactor 
complex as a function of neutrino energy are plotted for two distances to SADO, 50 km (blue) and 60 
km (red). The upper and lower boundary of each band corresponds, respectively, to sin 2 6\2 = 0.25 
and 0.35. Here Am^ = 7.9 x 10 -5 eV 2 and #13 = 0. The vertical dashed line marks the current 
KamLAND energy threshold E u = 3.4 MeV for reactor neutrino analysis, which corresponds to 
the maximum energy for geo-neutrinos. 

(A) SADO with single reactor complex, which will be denoted as SADO S i n gi e . 

(B) SADO with multiple reactor complex, which will be denoted as SADO mu i t i. 

The set up (A) may be an excellent approximation for many reactor sites in the world 
including Angra reactor in Brazil and Daya Bay in China [13]. On the other hand, the 
set up (B) with multiple reactor complex in addition to the nearest one, which could act 
as "background" , is inevitable if we think about experiments in Japan, for example, where 
many nuclear power plants (NPPs) are located within ~ 100 km. 

For definiteness, we take a particular site for the set up (B) in this paper; SADO at Mt. 
Komagatake, Niigata prefecture, Japan, which is about 54 km away from the Kashiwazaki- 
Kariwa nuclear reactor complex with the maximal thermal power of 24.3 GW t h, the largest 
NPP in the world. 2 Taking the particular setting is necessary to estimate the "background 



2 The site at ~ 1000 m below the peak of Mt. Komagatake is tentatively selected because of the appropriate 
distance from the Kashiwazaki-Kariwa NPP, a sufficient overburden, and possibility of relatively short 
access of about 3-5 km from a public road by digging a tunnel. Nonetheless, the site is not meant to 
be unique and there exist many other potential sites in the strip of distance 50-70 km from the reactor 
because the region is quite mountainous. 
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effects" caused by the other reactors than the nearest one. For the former case (A) one can 
use any reactor site, but for purpose of comparison, we just switch off contributions from the 
other 15 reactors except for the Kashiwazaki-Kariwa NPP to use the common normalization 
factor. In Table HVI in Appendix El we present the names of the reactor NPPs, their current 
(near future) thermal power, their distances from KamLAND and from SADO as well as 
their relative contributions to the neutrino flux at these sites. 

To make our analysis useful for experiments at any other reactor sites we present our 
plots, except for FigJ21 in units of GW tn -kt-yr where "kt" is an abbreviation of kton. Notice 
that the number in the unit refers only to the reactor power of the principal reactor, the 
Kashiwazaki-Kariwa NPP in this paper, and does not include the one from other 15 reactors. 
In this work, we ignore all the Japanese research reactors as well as any other reactors outside 
Japan. In fact, their contribution to the u e flux at KamLAND is calculated to be 4.5% 0, 
while the corresponding value at SADO is significantly smaller, about 1.1% due to the 
much larger flux from the closest reactor. We do not take into account these contributions 
explicitly in our calculation assuming that they can be subtracted from data. We can argue, 
quite reasonably, that the uncertainty of estimation of such contribution is at most 10-20%, 
and the additional systematic error due to the subtraction is negligibly small. 



C. KamLAND vs. SADO 

One of the crucial questions for us is "to what extent is the principle of tuning the 
experiment to the oscillation maximum effective in improving the accuracy of determination 
of 0i2?" The answer to this question will also tell us if SADO can supersede KamLAND, 
and if yes, to what extent. 

To answer this question, we present in Fig. El a comparison between the sensitivities to 
sin 2 6*i2 attainable by SADO sing i e , SADO mu i ti located at 54 km, and KamLAND. This plot 
shows how accurately we can reproduce the input mixing parameters, sin 2 9 12 = 0.29 and 
Anijj = 7.9 x 10~ 5 eV 2 , which will be used throughout this paper. 

For the purpose of comparison we have used for both, SADO and KamLAND, the same 
systematic error of 4% and the common energy threshold of prompt = 0.9 MeV. We deter- 
mined the exposure time of the each detector so that they receive (approximately) the same 
number of events without oscillation. By equalizing the number of events in this compar- 
ison, we aim to reveal how efficient the principle of tuning to the oscillation maximum for 
determining 8u is. For SADO mu i t i and KamLAND we have included all 16 NPPs, while for 
SADO S i n gi e we have considered only contributions from the Kashiwazaki-Kariwa NPP. (The 
distances to the 16 NPPs from SADO and KamLAND are given in Table ITV1 ) 

In comparing KamLAND with SADO in this subsection, we have ignored the geo-neutrino 
contributions. Therefore, the comparison between SADO mu i t i and KamLAND is purely a 
comparison of the two different detector locations. The possible impact of geo-neutrinos on 
SADO will be described in Sec. IIII1 and details of the analyzes in Sees. IIVI and IVl 

We discover in Fig. El the enormous power of SADO; the sensitivity attainable by SADO 
is better by a factor of ~ 2.5 than that of KamLAND for about three times longer operation 
in kt-yr. We also notice that the sensitivities achievable by SADO S i ng i e and SADO mu i t i are 
comparable; the additional contribution from the other 15 reactors does not affect so much 
the #i2 determination. More quantitatively, we recognize from Fig. El that at 90% CL, the 
sensitivity attainable by SADO sing i ejmu iti corresponds to about 7% (5%) error in sin 2 #i 2 for 
0.5 kt-yr (1.5 kt-yr) measurement. 
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FIG. 2: Comparison between the sensitivities at 90% CL for 2 DOF expected by KamLAND and 
SADO for runs with (about) equal number of events without oscillation. Neutrinos from all the 
16 NPPs are taken into account for KamLAND whereas for SADO s i ng i e only Kashiwazaki-Kariwa 
NPP contribution was computed. The input values of mixing parameters are sin 2 812 = 0.29 and 
Amli = 7.9 x 10~ 5 eV 2 . The numbers of events for 0.5 kt-yr measurement at SADO s i ng i e and for 
1.4 kt-yr measurement at KamLAND will be about 1700. The systematic errors are assumed to be 
4% in both detectors. Events with prompt energy (-E pr0 mpt = E u - 0.8 MeV) greater than 0.9 MeV 
are used for both detectors. For a comparison, we also show the sensitivity, at 90% CL, expected 
for SADO mu iti if the contributions of all 16 NPPs were included. 



One may want to ask a further question. Namely, how long does KamLAND take to reach 
the same sensitivity as of SADO? (In fact, this is the question raised by one of the referees.) 
The answer to the question is in fact a surprising one: To achieve the same sensitivity of 
SADO mu iti presented in the upper panel of Fig.|21 KamLAND would need more than 100 kt-yr 
exposure. Whereas for the cases SADO S i ng i e shown in the upper panel and SADO S i ng i e ,muiti 
shown in the lower panel, KamLAND would not achieve the same sensitivities even in the 
limit of infinite statistics. Therefore, tuning the baseline distance is essential to achieve the 
ultimate sensitivity of 2% level for sin 2 9u determination. 

We should mention that, for the lower value of Am^ = 7.2 x 10~ 5 eV 2 , preferred by the 
data |5j before KamLAND 2004 result, the optimal distance would have to be re-scaled to 
about 70 km (see Eq.Q), resulting in a corresponding increase of a factor of two in running 
time in SADO to match the same sensitivity at higher Am^ achieved for 54 km. 
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D. sin 2 6>i2 versus sin 2 2#i 2 



Though it is slightly off the main line of the discussion let us make some comments on 
usage of the variable in the plots we present in this paper. There have been various choices 
of variables for displaying allowed regions in # 12 - Am^ plane. They include, tan 2 6* 12 , 
sin 2 26^2, and sin 2 6 l 12 . We feel it important, if possible, to have a "standard choice" of the 
mixing variable to make comparison between different analysis easier. We choose sin 2 # 12 
throughout this paper for reasons which are described below. 

First of all, we now know that 9 12 is in the first octant. Therefore, we have no more 
reasons to continue to use tan 2 6*i 2 . Then, the question is which of sin 2 26*i 2 , or sin 2 9 12 is to 
be used. This is a difficult choice because both of them can be supported for good reasons 
from physics point of view. The solar neutrinos experience neutrino flavor transformation by 
different mechanisms at low and high energies. At low energies, it is essentially the vacuum 
oscillation effect that is dominant. At high energies the adiabatic MSW conversion effect 
takes over where the u e survival probability is approximately given by sin 2 6* 12 . Therefore, 
the natural variable at low and high energies are sin 2 26 l 12 , and sin 2 # 12 , respectively, and 
they both have right to be chosen. 

We use sin 2 9 12 in this paper because of two reasons. One is that sin 2 9\2 allows a direct 
physical interpretation as the probability of finding v e in the second mass eigenstate, as em- 
phasized in [27| . The other one is that #i 2 is relatively large, though not quite maximal. For 
large mixing angle 9, sin 2 29 is not a convenient variable because sin 2 29 becomes insensitive 
to change in 9. It can be simply understood by computing the Jacobian 

5 (sin 2 29) d sin 2 29 
d(sm 9) asm 

where ^(sin 2 ^) and 5 (sin 2 29) denote the error in sin 2 9 and sin 2 29, respectively. It means 
that a large error obtained for sin 2 9 can be translated into a much more suppressed error for 
sin 2 29 because of the small Jacobian at large 9 ~ tt/4. Of course, the problem of disparity 
between errors estimated for these two variables in our case is much milder than that of 
atmospheric angle 923 [H 



III. REQUIREMENTS AND POSSIBLE OBSTACLES FOR THE PRECISION 
MEASUREMENT OF 12 

In order to achieve the optimal sensitivity promised in the previous sections the following 
two requirements have to be met. They are: (i) experimental systematic error of 4% level, 
(ii) baseline distance of 50-70 km. We will discuss the requirement (i) and other related issues 
in this section. We calculate the optimal distance after setting up our analysis procedure in 
the next section. 



A. Experimental systematic error 

The key ingredients in the accurate measurement of #i 2 in reactor experiments is the 
systematic error in the measurement. In this paper, we assume that there exists a dedicated 
front detector. It would be ideal if there already exists a near-far detector complex for reactor 



8 



measurement of #13 at less than 1-2 km from the principal reactor for SADO. Presence of 
the front detector is crucial to suppress the systematic errors for accurate measurement of 
9±2 because the various errors such as absolute flux and cross sections largely cancel between 
the front detector and SADO, in complete parallelism with the reactor 13 measurement. 
See e.g., (28l . 29l | for details on how and which errors can be canceled. 



Since SADO cannot have structure exactly identical with the front detector, the system- 
atic error may not be as small as ~ 1%, the number expected to be reachable in reactor # 13 
measurement k|. We assume in our analysis, based on [30], that the systematic error can 
be as small as 4% at SADO. 3 This systematic error is supposed to include background from 
cosmic-ray events, assuming that the overburden is enough. Our claim that 4% systematic 
error may be in reach can be regarded as reasonable by recalling that the CHOOZ exper- 
iment already achieved the goal of the systematic error (1.5% for detection efficiency 
error) apart from the one due to flux normalization. Therefore, we think that the 4% is a 
conservative estimate and a smaller systematic error of ~2% may well be thinkable. 

One of the consequence for the requirement of less than 4% systematic error is that we 
cannot rely on the event cut of the energy spectrum of prompt energy -E pr0 mpt > 2.6 MeV to 
reject geophysical neutrino contamination, as done by the KamLAND group in their reactor 
neutrino analysis (H E| . Because of the uncertainty in the energy measurement it produces 
an additional ~ 2% systematic error which makes over-all 4% systematic error difficult to 
achieve. Therefore, we must take all the events with prompt energy greater than 0.9 MeV. 
But it means that we have to deal with geophysical neutrinos in our analysis. It will be one 
of the most important issues in our sensitivity analysis in this paper. 



B. Geophysical neutrinos; a brief review 

Geo-neutrinos comes into play in our sensitivity analysis of 6*i 2 because U (Th) decay series 
processes through six (four) /3-decays, producing 6 (4) geo-neutrinos with energy E v < 3.4 
MeV, or E pmmpt < 2.4 MeV. The threshold of = 1.7 MeV at SADO then implies that 
geo-neutrino events from U and Th decays can contaminate reactor neutrino events in our 
analysis. 

Then, what are geo-neutrinos and do we know their properties? In the rest of this 
subsection, we briefly review what we know about geophysical neutrinos. The radioactive 
decay chains of radionuclides such as 40 K, 238 U and 232 Th inside the Earth produce not 
only heat but also electron antineutrinos. The fact that the radioactive heat flux and 
the geophysical neutrino flux are tightly linked was first pointed out by Eder j3jj in the 
sixties. Since then, the geo-neutrinos have been the subject of continuing interests j32|. 
Unfortunately, a quantitative accounting of the radioactive heat flux requires a detailed 



3 In the first and the second versions of this paper posted on the arXiv (hep-ph/0407326), we have used 
Amji = 8.2 x 10~ 5 eV 2 based on the earlier version of Q, which is 3.8% larger than the value used in the 
third and this published (fourth) version. The readers who want to examine how different are the results 
with slightly different values of Am 2 ^ can look into the second version of the paper. More importantly, 
we assumed systematic error of 2% (4%) in the first (second, third and fourth) version. Therefore, the 
readers who want to do a detailed comparison between the cases of systematic error of 4% and 2% are 
advised to look into the first version of the paper. The reason why the sensitivities to 612 are essentially 
unchanged is that the measurement is not systematics dominated even at 60 GWth-kt-yr. 
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knowledge of the abundance distributions of these long-lived radioactive species inside the 
Earth. Our current knowledge of these abundance distributions is incomplete because of 
many reasons, for example, direct sampling is possible only at or near the surface and most 
of the Earth's surface is inaccessible. Consequently, the radiogenic heat flux and therefore 
the geo-neutrino flux is largely model-dependent and its precise magnitude is unknown. 

The presently estimated value for the global heat flux from the earth is 40 TW. It gives 
an upper bound on geo-neutrino flux from the interior of the earth. Heating from known 
radioactive sources in the surface layers (corresponding to 1/300 of the Earth's total mass) 
is estimated to be already about 20 TW. Most geophysical models predict that the concen- 
trations of radioactive sources will decrease rapidly with depth, but even small variations in 
these concentrations, for such a big mass as that of our planet, can greatly affect heat and 
geo-neutrino production. 

We will treat geo-neutrinos flux as an additional parameter to be fit in our sensitivity 
analysis. To take proper care of its uncertainty we examine the two extreme cases of no 
geo-neutrino flux as an input and the maximum flux corresponding to global heat flux of 
40 TW (the Fully Radiogenic Model). Fortunately, as we will see in the next two sections, 
the sensitivity on 6*12 to be achieved by SADO is not disturbed in an essential way by the 
presence of geo-neutrinos provided that we choose appropriate baseline distances. 



C. Uncertainty due to 



'13 



Uncertainty due to the unknown value of #13 can be another source of the systematic 
error. It was shown in |33( that it can be regarded as an effective uncertainty in the flux 
normalization. Let us briefly review their argument. Since the vacuum neutrino oscillation 
is a good approximation for the reactor 9 12 experiments, we can write down the electron 
antineutrino survival probability in a simple approximate form 



P(v P 



sin #i3 + cos 9 



13 



-I f ■ 2nn ■ 2 ( Am^L,- 

1 — Ji sm sin 



V 4E„ 



(4) 



where E u is the neutrino energy, Li is the distance to the detector from reactor i, and 
fi is the fraction of the total neutrino flux which comes from reactor %. Thus, nonzero 
#13 effectively acts as an uncertainty in the flux normalization of less than 8%, giving the 
CHOOZ constraint sin 2 29 13 < 0.15 at 90% CL Q. 

This uncertainty in the flux normalization cannot be canceled by measurement at the 
front detector. Therefore, #13 has to be determined by an independent measurement. But 
the experimental determination of # 13 should always come with errors, 5(sin 2 2# 13 ). It will 
produce an effective uncertainty in the probability by the amount 



S(cos 4 9 13 ) 



d cos 4 9 



13 



d(sin 2 29 



13 



S(sm 2 20i 3 ) 



1 + tan 2 9 



13 



5(sin 2 20 13 ), 



(5) 



for a given experimental error £(sin 2 29 is) of measurement. If we take the estimation in 
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for reactor (9 13 experiment, 5(sin 2 2# 13 ) ~ 0.012 (0.02 at 90% CL) almost independently 
of the true values of sin 2 2# 13 . Therefore, we expect an uncertainty of ~ 0.6% level in the 
probability and hence in sin 2 2# 12 . It may be translated, by using (jSJ), into an uncertainty 
on sin 2 9 12 of about 0.4% at the current best fit of #12. 
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Thus, the error of #13, once it is measured by reactor experiments, adds only ~ 0.4% 
to the uncertainty on sin 2 612. Though it does not pose any serious problem, more precise 
measurement of #13 either by LBL high-intensity superbeam experiments with both neutrino 
and antineutrino modes 0, 17, 3], or by high statistics reactor experiments or even 



by neutrino factory |35( is highly welcome. 

We assume in the rest of this paper vanishing # 13 for simplicity. Assuming 13 is deter- 
mined with error much smaller than 1% the effect of nonzero #13 effectively acts as smaller 
v e flux by a factor of cos 4 #13. For example, if #13 is on the edge of the above CHOOZ bound, 
cos 4 #13 = 0.92, we need 9% longer running time to obtain the same sensitivities as presented 
in the later sections. 



IV. ANALYSIS PROCEDURE AND ESTIMATION OF THE OPTIMAL DIS- 
TANCES 

In this section, we first make a rough description of our analysis procedure, and then 
estimate the optimal distances for reactor #12 experiments. For more details of the analysis 
procedure, we refer the readers to Appendix IH1 



A. Analysis procedure in brief 

We follow our previous papers [3^, |3tJ in calculating the number of neutrino events 
expected from reactors as well as from Th and U radioactive decays in the Earth at a liquid 
scintillator detector, SADO. We will give our results in terms of GW t h ■ kt ■ yr. Notice that 
10 GW t h ■ kt ■ yr is approximately equivalent to 1 year exposure of 0.5 kt detector at 54 km 
away from Kashiwazaki-Kariwa NPP operating at 80% efficiency 

The total number of events from a single NPP for an exposure of 1 GW t h ■ kt ■ yr , for 
the used threshold energy is estimated to be, 

-2 

[GW th -kt-yr]-\ (6) 

where the upper (lower) value corresponds to no oscillation (oscillation at the input values 
used in this work). GW t h is defined to be the thermal power actually generated, which 
should not be confused with the maximal thermal power of a given NPP. 

The total number of events from geo-neutrinos expected from the Fully Radiogenic 
Model |32j at Mt. Komagatake is 

AC al (£ prompt > 0.9 MeV) = 49 [kt • yr ]~\ (7) 

Since we currently do not know the U and Th geo-neutrino fluxes we consider two extreme 
cases: (i) zero geo-neutrino input and (ii) geo-neutrino input flux from U and Th computed 
assuming that the radiogenic production accounts for the total Earth heat flow of 40 TW 
(input from the Fully Radiogenic Model). They represent two extremes of the input geo- 
neutrino flux between zero and 3xl0 7 cm -2 s" 1 at the position of the detector, and one 
would expect the real situation to be somewhere between these two numbers. Note that 
even if we assume zero geo-neutrino input, we allow non-zero geo-neutrino flux as an output 
in our fit which can confuse reactor neutrinos. Since the ratio of U to Th contributions is 



N™(E m t > 0.9 MeV) 



141 

46 



54 km 
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rather well defined and common in varying models, we assume in our analysis that their 
relative contribution is held fixed while we treat an over-all absolute normalization as a free 
parameter. The Th contribution is assumed to be 83% of the total geo- neutrino flux [32]. 
This procedure should be (and will be) tested by future KamLAND data which includes 
-E'prompt < 2.6 MeV. Of course, having the additional free parameter in general makes the 
sensitivity on # 12 worse. 

We use the following definition of the x 2 function as 

X (sin 0i2,Am 21 ,0u) = V- g ' ( 8 ) 

^ a t 

i=i ' 

where cr, = \J N^ xp + (0.04 N°* p ) 2 is the statistical plus systematic uncertainty in the number 
of events in the i-th bin; Nf is the theoretical expected number of events as functions of 
mixing parameters as well as <j)\j, the total geo- neutrino flux from U, to be fitted, calculated 
as explained in the Appendix lB| A^ exp = A^ heo (sin 2 6*12 (input), Am^ (input), 0tj (input)) is 
the simulated expected number of events at the detector for the input values of the mixing 
parameters as well as geo-neutrino flux, assumed in this paper, i.e., sin 2 6*12 (input) = 0.29 
and Am^ (input) = 7.9 x 10~ 5 eV 2 and the geo-neutrino input, 0u(i n put), will be either 
zero or equal to what is expected if the radiogenic contribution to the terrestrial heat is to 
be explained by the Fully Radiogenic Model. The crustal contribution was estimated as in 
Ref. |3j| but for a detector located at Mt. Komagatake in Niigata prefecture, Japan. 4 

Using the x 2 function we compute the region in the (sin 2 #12, Am^) space allowed by 
SADO spectrum data at 68.27%, 90%, 95%, 99% and 99.73% CL for a given GW th -kt-yr 
exposure. 



B. Optimal baseline distance 

We have discussed in Sec. II that it is the key for a precise measurement of 9 12 measure- 
ment to choose the baseline distance which corresponds to the oscillation maximum. To 
seek ultimate sensitivity on #12, however, we must first elaborate this point. 

We show, in the upper and the lower panels in Fig. the accuracy of determination 
of sin 2 6*12 and Am^, respectively, as a function of L expected by 10 GWth-kt-yr exposure 
computed by the minimization of the x 2 function that is defined in Sec. III. The results 
for two different energy thresholds, 0.9 and 2.6 MeV, are presented for the purpose of 
comparison. For the case -E pr0 mpt > 0.9 MeV, we show two curves for sensitivity with and 
without geo-neutrino background contribution. We note that both Am^ and 9 12 are left free 
and marginalized in the analysis to obtain the sensitivity of sin 2 # 12 and Am 2 1; as presented 
respectively in the upper and the lower panels in Fig. |H1 

We see that when the analysis threshold is higher the best place to put a detector, in 
order to achieve greater sensitivity to #12, is at around L ~ Lou = 63 km, in agreement 
with our rough expectation. On the other hand, when the threshold is lower, it is preferable 
to choose values of L 10-25% smaller than Lqm 5 - 



4 Be careful about the fact that there are at least 13 Mt. Komagatake in Japan. 

5 This fact was recognized in a careful study of sensitivity in reactor 13 measurement |lj2|2j|. But, the 
shorter baseline L by 10-25% does not make a great difference in that case because L is short, ~ 1 km. 
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Distance (km) 

FIG. 3: Achievable precision of the determination of the oscillation parameters at 68.27% CL for 1 
DOF as a function of the detector distance from the Kashiwazaki-Kariwa NPP for 10 GWth - kt-yr 
exposure and energy thresholds of 0.9 and 2.6 MeV. The horizontal dashed lines indicate the input 
values of the mixing parameters. The geo- neutrino contribution was calculated as predicted by the 
Fully Radiogenic Model. For the lower energy threshold the x 2 was minimized allowing the total 
geo-neutrino flux to vary from to 3xl0 7 cm~ 2 s^ 1 . We have used the input values of the mixing 
parameters corresponding to sin 2 9\2 = 0.29 and Am^ = 7.9 x 10 -5 eV 2 . The vertical dashed lines 
indicate the optimal range for SADO detector. 

It is remarkable that if there were no geo-neutrinos the baseline distance as short as 40 
km could be suitable. Even when we take them into account, the optimal distance ranges 
between 50 and 70 km which is still a relatively wide range. We will see below (Sec. |V} that 
even at 70 km, the longest end of the distance range, the sensitivity is worse only by about 
15% or so than SADO at 54 km. This point is very important and encouraging for people 
who try to design reactor 6*12 experiments because a wide band of radii between 50 and 70 
km around a particular reactor NPP can serve as a potential location for the experiment. 

The reason why the effect of geo-neutrinos is larger at shorter distances is as follows. As 
L gets shorter the region where oscillation effect is sizable moves to lower energies. Then, 
for L < 50 km it starts to invade the region of E u < 3.4 MeV where geo-neutrinos live. 

An interesting behavior of the error on sin 2 9%2 for the 2.6 MeV threshold, rapidly increas- 
ing with L from about 80 to 150 km and decreasing again from 150 to about 170 km, can be 



But in reactor 812 experiments, however, the difference is significant because the choice of L = 50 km or 
L = 70 km implies a completely different site for SADO. 
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explained as follows. The maximal sensitivity, reached at about 60 km, corresponds to the 
first oscillation maximum at the peak of the event distribution. As L increases this first os- 
cillation maximum moves to higher energies, finally going out of the reactor neutrino energy 
spectrum. When the second oscillation maximum passes the 2.6 MeV threshold, at about 
150 km, it improves the sin 2 # 12 determination causing the partial recovery of precision. 

In the lower panel of Fig. |H] we appreciate the weak dependence on L of the accuracy in 
the determination of Am^ for L > 50 km. We will see in Sec. IV11I fTable llllJ) that SADO 
does improve the sensitivity of Am 21 determination over the KamLAND's though not to 
the extent that occurs for Qu- 
it is interesting to clarify the difference between our case and that of the experimental set 
up of the next generation LBL accelerator neutrino experiments 0, 17, la] . In particular, 
in the JPARC-SK experiment, the accuracy of determination of sin 2 2# 23 is expected to 
reach down to ~ 1% at 90% CL [16]. In this experiment, where L is already fixed, the 
best sensitivity for the angle determination is obtained when the experiment is tuned to 
the energy which corresponds to the first oscillation maximum. The spectrum of an off- 
axis beam is nearly monochromatic. On the other hand, in the reactor experiment we are 
considering, the energy spectrum is given and is quite broad and hence nearly maximal 
sensitivity to the mixing angle can be achieved in a range of values of L. 



V. HOW ACCURATE IS THE REACTOR MEASUREMENT OF 12 ? 

In this section, we give a quantitative estimate of the sensitivity which can be achievable 
in a dedicated reactor measurement of 6*i 2 . 



A. Analysis of sensitivity on 6>i 2 : SADO at 54 km 

We show in Fig. HJhow precisely the mixing parameters can be determined by SADO S i ng i e , 
assuming that only the Kashiwazaki-Kariwa NPP will be contributing to the neutrino flux 
at the detector site, for exposures of 10 GWth-kt-yr (upper panels), 20 GWth-kt-yr (middle 
panels) and 60 GWth-kt-yr (lower panels). In the left panels, geo-neutrinos were not consid- 
ered, in the middle panels, we have set the geo-neutrino flux input to be zero but allowed 
the flux to be fitted to vary freely in the \ 2 fit- In the right panels, we have calculated 
the geo-neutrino input according to the Fully Radiogenic Model and allowed the flux to be 
fitted to vary freely in the \ 2 fit- 

In Fig. we provide the same information as in Fig. |U but for SADO mu i ti which includes 
the contributions of all the 16 reactor NPPs in Japan. We note that by taking into account 
the contributions from the other 15 reactor NPPs, the accuracy of determination of sin 2 9\2 
becomes slightly worse. The error on Am^ also slightly increases. 

The inclusion of geo-neutrinos does not influence very much the sensitivity to sin 2 9 12 as 
can be seen in Figs. 0]and|3J This is because we have chosen a distance such that the position 
of the dip is far enough from the geo-neutrino energy region. This can be visualized in Fig. 
|BJ where we present the normalized expected energy spectra at SADO sing i e (left panel) and 
SADO mu iti (right panel) for exposures of 10 and 60 GW tn -kt-yr. The impact of the inclusion 
of the 15 NPPs is to modify somewhat the shape of the spectrum. 

On the other hand, the sensitivity to Am 21 is slightly more affected by the inclusion of 
the geo-neutrino background. This is because in Am 21 determination, it is important to 
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FIG. 4: Precision of the determination of the oscillation parameters by SADO s i ng i e , a detector 
at 54 km from the Kashiwazaki-Kariwa NPP, for exposures of 10 GW t h - kt-yr (upper panels), 20 
GWth-kt-yr (middle panels) and 60 GW t h-kt-yr (lower panels). Left and middle panels have no 
geo-neutrino input, but in the middle panels geo-neutrinos were however taken into account in the 
fit. Right panels were calculated for geo-neutrino input according to the Fully Radiogenic Model. 
The confidence level regions were computed for 2 DOF. 
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FIG. 5: Same as Fig. |1] but for SADO mu i t i with all 16 reactor NPPs included. 



observe the shape of the spectrum in the entire energy range making it more sensitive to 
the geo-neutrino contributions in the lower energy bins. At this point it is also instructive 
to look at Fig. [3 where we show the A% 2 = x 2 behavior as a function of sin 2 9 12 and Am 
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for 10, 20 and 60 GW t h-kt-yr exposure, all 16 NPPs considered. 

We summarize in Table [Q the accuracy of determination of the mixing angle 9 12 which 
can be achieved by SADO S i ng i c and SADO mu i t i for these three exposures. We use the prompt 
energy threshold of 0.9 MeV and take into account the background effect of geo-neutrinos 
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FIG. 6: Spectral distortion expected as a function of the e + prompt energy for SADO s i ng i e (left 
panel) and SADO mu i t i (right panel). Each ratio is represented as a point with the corresponding 
error bar for a given exposure (10 and 60 GW t h - kt-yr). The number of oscillation events at each 



bin was calculated for the input values sin 612 



0.29 and Am^ 



7.9 x 10" 5 eV 2 . The dashed 



line indicates the maximal energy for geo-neutrinos. Geo-neutrino contribution from the Fully 
Radiogenic Model was included. 



TABLE I: The expected sensitivity to the solar mixing angle at SADO s i ng i e (SADO mu i t i) at 1 a 
(68.27%) and 3 a (99.73%) CL for 1 DOF obtained for SADO at Mt. Komagatake (54 km, the 
upper rows) and Sado gold mine (71 km, the lower rows). They are calculated with the prompt 
energy threshold of 0.9 MeV and with the background effect of geo-neutrinos from Fully Radiogenic 
Model. The numbers for SADO mu i t i are computed for both detectors by using position-dependent 
contributions from all the 16 NPPs. 



exposure 


sin 2 6>i2 [68.27% CL] sin 2 12 [99.73% CL] 


54 km 


10 GWth-kt-yr 


n 9Q +0.014 (( , 9Q +0.015n n 9Q +0.040 / n oq +0.042n 
U - 2y -0.013 l U - 2y -0.014J U - 2y -0.031 l U - 2y -0.033J 


20 GWth-kt-yr 


n 9Q +0.010 (n. 90+O.Oin n 90+O.O28 / n 9Q+0.03CK 
U.29_ 0010 (0.29_ 0011 ) 0.29_ 0023 (U-29_ 025 ) 


60 GWth-kt-yr 


n 90+O.OO6 / n 2 q+0.0(m 90+O.OI7 / 2 n+0.019x 
u - zy -0.006 V u -^ y -0.007^ u - zy -0.015 V u - zy -0.017J 


71 km 


10 GWth-kt-yr 


q 90+O.OI6 / 2 q+0.019N 9O+0.045 ((] 2 q+0.052n 
u - zy -0.015 V u -^ y -0.017^ u - zy -0.037 V u - zy -0.040J 


20 GWth-kt-yr 


q 9O+0.012 / n 2 o+0.013n 9O+0.031 / 9 q +0.036N 
u - zy -0.011 V u -^ y -0.012^ u - zy -0.027 V u - zy -0.030J 


60 GWth-kt-yr 


9O+0.007 fn 2O +0.008n Q 2q +0.018 ((] 2q +0.022x 
u - zy -0.007 V u -^ y -0.008^ u - zy -0.016 V u - zy -0.020J 
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FIG. 7: Ax 2 as a function of the oscillation parameters around the minimum for 10, 20 and 
60 GW t h-kt-yr exposure. All 16 reactor NPPs were included in the calculation and the prompt 
energy threshold used was 0.9 MeV. Geo-neutrino contribution from the Fully Radiogenic Model 
was included. 

expected by the Fully Radiogenic Model, the most conservative estimate. We have investi- 
gated the dependence of this accuracy on the input value of Am^ by varying it in the range 
Amjj = (7 — 9) x 10~ 5 eV 2 . We have discovered the accuracy becomes somewhat worse, at 
most by 30%, if the input value is smaller than the current best fit one due to geo-neutrinos, 
whereas for higher values there is practically no change. 

In Table HI we also present for comparison the results of our analysis for a detector placed 
at Sado island, about 71 km away from the Kashiwazaki-Kariwa NPP. As one can anticipate 
from Fig.|3]the sensitivity at 71 km is worse than that of SADO at 54 km but only by 15-20%. 
Thus, SADO at the Sado island is still a valid option for the reactor 9 12 experiment. 

VI. STABILITY OF THE RESULTS AGAINST CHANGES IN STATISTICAL 
PROCEDURE 

We have shown in the previous sections that the sensitivity for 9 12 attainable by SADO 
with a modest requirement of ~ 4% systematic error is extremely good, which is competitive 
to all the methods so far proposed. In this section we confirm the stability of the results 
by using different definitions of x 2 ■ While we believe that our statistical method and the 
treatment of errors based on x 2 defined in Eq. (JSJ) is reasonable and a standard one, it is 
nice if we can explicitly verify that our results are stable against changes in the statistical 
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procedure. 

For this purpose, we examine the following two possible choices of x 2 ■ The first one, 
which is frequently used in various statistical analyzes of data, reads 



2 A(ATf p -«iV7 he °) 2 (a-1) 2 

xLr = nun N c xp + a2 ' , (9) 

i=1 i syst 

where iV 4 exp and j\T theo are defined in the same way as in our x 2 definition in ©. The x 2 
was used, for example, in the KamLAND analysis in Ref. [38|. We take, as in our analysis, 
cr syst = 4%. As we will see below, this choice of x 2 leads to even higher sensitivity of 6*12 
measurement at SADO. It is not surprising to obtain such result because the systematic 
uncertainty affects all bins simultaneously for x 2 in © ; whereas it can fluctuate bin- by-bin 
with our definition of x 2 m Eq. (JHJ). 

Our second choice of x 2 is the one used for the analysis in reactor #13 measurement (28^ . 
In such experiments we compare yields at identical near and far detectors and a large portion 
of the systematic errors cancels out in such a setting. Therefore, the statistical treatment 
involves the errors which have different characteristics, (un-) correlated between near and far 
detectors, and between bins. Following jH, , we consider four types of systematic errors: 
c"db, c"Db, 0"dB, and cr dh . The subscript D (d) represents the fact that the error is correlated 
(uncorrelated) between detectors. The subscript B (b) represents that the error is correlated 
(uncorrelated) among bins. The definition of x 2 is 



xli = min V 

as ^— ~ ' 

a=f,n 
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„2 I „,2 



(Ageo _ {1 + ai + aa + a)N »*y 



a 2 



(10) 



where N^ eo represents the theoretical number of events at the near (a = n) or far (a = /) 
detector within the z-th bin. Again, N^ p is defined as the number of signal event calculated 
with the best-fit parameters of the "experimental data" . See [29] for more details. 

To simulate the reactor #13 measurement with sensitivity up to ~ 1% which is, very 
roughly speaking, equivalent to sensitivity up to sin 2 26 ) i 3 = 0.01 for a long enough exposure 
(systematics dominated measurement), the following numbers have been taken for these 
errors ctdb = o"Db = 2.5% and o^b = o"db = 0.5%. To simulate the systematic error of ~ 
4% at SADO we tentatively multiply by 4 all these errors. That is, we take o"db = o"Db = 10% 
and <JdB = o"db = 2% in our analysis for SADO. We also examine the optimistic case of 
systematic error of ~ 2% at SADO, considering the possibility that the systematic error can 
be improved. The errors for this case are taken as 1/2 of the 4% case. We have confirmed in 
the case of #13 measurement that the error twice (four times) larger than that of [2£| leads 
to the sensitivity approximately equal to sin 2 26*i3 = 0.02 (0.04) for long enough exposure. 

Since the front detector and SADO will be different in volume, assuming equal errors 
for both detector is nothing more than a crude approximation. But, we feel that it gives a 
reasonable framework to check the stability of our statistical treatment. 

In Fig. [S| we show the result of our analysis with three different definitions of x 2 1 © and 
(fTUJ) defined above, and the one given in JHJ) that is adopted in our analysis. We observe that 
our result obtained in the previous section by using the x 2 in © lies between the results 
of the analysis with two different x 2 1 in © an d (|10|). It is expected that the result with 
X 2 in © gives tighter constraint on the oscillation parameters because the way how the 
systematic error is treated only allows fluctuation of the absolute normalization, not bin-by- 
bin independent fluctuations, and it is harder to mimic spectral shape distortion. We do not 
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FIG. 8: Regions allowed at 99.73% CL for 60 GW t h-kt-yr exposure, determined by using three 
different \ 2 definitions: the definition used throughout this paper (Eq. (JHJ)), Xcorr defined in Eq. @ 
and defined in Eq. (|1U|) . The contours with \ 2 defined above are, in order, the middle (black 
line), the thinnest (red) and the thickest (blue) ones in widths. 



have any good physics intuition on what would be the result for the choice of x 2 ; in (|10|b 
but it turns out that the sensitivity is a little worse than our estimate with (JHJ). Nonetheless, 
the difference between the results obtained by using three different definitions of \ 2 is not 
significant. Hence, we conclude that our estimate of the sensitivity of Qyi determination at 
SADO is stable under change of statistical treatment. 



VII. COMPARISON OF SENSITIVITIES OF REACTOR 12 MEASUREMENT 
WITH OTHER METHODS 



In this section, we compare the sensitivity calculated in Sec.|3for a reactor measurement 
of #12 at SADO with the ones expected by other methods, in particular the solar neutrino 
experiments. We assume CPT symmetry in the discussions in this subsection. 

To date the most accurate determination of Qyi can be accomplished by combinin g al l 
the solar neutrino experiments and KamLAND data 0]. According to Table 2 of Ref. |38|. 
the current solar neutrino data together with KamLAND result can be used to determine 
sin 2 #12 with about ~ 8% precision at 68.27% CL for 1 DOF. From the Table 3 of Ref. Q, 



one expects that this precision will not be significantly reduced even if one combines with 
future KamLAND data corresponding to 3 years of exposure. 
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TABLE II: Comparisons of fractional errors of the experimentally determined mixing angle, 
6s\ 2 /s\ 2 = <5( srn2 $i2)/ sin 2 #12, by current and future solar neutrino experiments and KamLAND 
(KL), obtained from Tables 3 and 8 of Ref. 141 ] . versus that by SADO s i ng i e (SADO mu i t i) obtained 
in this work at 68.27% and 99.73% CL for 1 DOF. 



Experiments 


5s\ 2 /s 2 12 at 68.27% CL 


8s{ 2 /s\ 2 at 99.73% CL 


Solar+ KL (present) 


8 % 


26 % 


Solar+ KL (3 yr) 


7% 


20 % 


Solar+ KL (3 yr) + pp (1%) 


4% 


11 % 


54 km 


SADO for 10 GWth-kt-yr 


4.6% (5.0%) 


12.2 % (12.9 %) 


SADO for 20 GW th -kt-yr 


3.4% (3.8%) 


8.8 % (9.5 %) 


SADO for 60 GW th -kt-yr 


2.1 % (2.4 %) 


5.5 % (6.2 %) 



It is expected that if future solar neutrino experiments can detect pp neutrinos selectively, 
they would give the best sensitivity to # 12 determination. Assuming the most optimistic 
error of 1% for pp neutrino measurement and by combining with the other solar neutrino 
experiments as well as with 3 years running of KamLAND, one expects, from Table 8 of 
Ref. 13 . that it is possible to measure sin 2 #i 2 with uncertainty of ~ 4% at 68.27% CL 



for 1 DOF. Despite the fact that this result is obtained with the previous best fit value of 
Am 21 (~ 7 x 10~ 5 eV 2 ) it appears to be safe to assume that it will remain almost unchanged 
because the final precision of Q\ 2 is essentially determined by the solar neutrino data. 

The sensitivities expected by a reactor neutrino experiment calculated in the previous 
subsection should be compared with the precision obtained by all solar + KamLAND experi- 
ments mentioned above. This is done in Table ILT1 from which one can see that SADO sm gie,muiti 
with exposure longer than 20 GWth-kt-yr, can determine #12 better than all other observa- 
tions combined. We want to mention that Gd-loaded Super-Kamiokande [3^ cannot com- 
pete with the SADO's sensitivity on 9u Even after combining with the solar data the 
error 15% on sin 2 9\ 2 at 99% CL for 3 years operation is larger than SADO's 12% at 99.97% 
CL for its 10 GW th -kt-yr (~ 0.5 kt-yr at SADO) exposure. It should be stressed that SADO 
alone can achieve such great precision without being combined with any other experiments, 
assuming that 6* 13 is measured with reasonable precision. Therefore, we conclude that a 
reactor measurement can supersede other methods for precise measurement of 6\ 2 if the 
experimental systematic error of 4% is realized. 

Although our main concern in this paper is a precise determination of #12 let us briefly 
discuss the sensitivity of SADO on Am 2l and its comparison with the one by other methods. 
In Table ITLT1 we present the expected sensitivity to Am 2l at SADO sing i (SADO mu i t i) at 1 
a (68.27%) and 3 o (99.73%) CL for 1 DOF, obtained with the prompt energy threshold 
of 0.9 MeV and with the background effect of geo-neutrinos from Fully Radiogenic Model. 
The present KamLAND data together with all the solar neutrino data already achieved 
14% error at 3a CL on Am 21 determination. As one can see from Table ITTT1 (and Fig. 
SADO S i n gie (SADO mu iti) will be able to reduce the error down to 4.3% (4.6%) at 3a CL for 
60 GWth-kt-yr operation. The error is smaller by a factor of 2 than that expected for the 
combined analysis of the future solar neutrino data and 3 years operation of KamLAND 
found in Fig. 6 of Ref.j3]. Nevertheless, it is still greater than the error of about 2.8% at 
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TABLE III: The expected sensitivity to the Amfi at SADO sing i e (SADO mu i ti ) at 1 a (68.27%) 
and 3 a (99.73%) CL for 1 DOF, obtained with the prompt energy threshold of 0.9 MeV and with 
the background effect of geo-neutrinos from Fully Radiogenic Model. 



54 km 


exposure 


AmJ, x 10 5 eV 2 [68.27% CL] 


Am 2 12 x 10 5 eV 2 [99.73% CL] 


10 GWth-kt-yr 


7.901^ (7.901858) 


7.9ot8:l (7.9018:?°) 


20 GWth-kt-yr 


7.9018;I1 (7.90l8lx) 


7.9018:^? (7.9oi8:^) 


60 GW th -kt-yr 


7-90±°it (7.9018:!!) 


7.90l8:| (7.901818) 



3a CL achievable by Gd- loaded Super-Kamiokande for its three years operation |40j. 

VIII. PHYSICS IMPLICATIONS OF ACCURATE MEASUREMENT OF 12 

There exist a variety of physics implications available when accurate measurement of 9 12 
is carried out. Here, we describe only a part of them. Of course, #12 is one of the fundamental 
parameters of particle physics and its precise determination is by itself clearly important, as 
already discussed in Introduction. It may be appropriate to add a remark in this context. 
Namely, the determination of sin 2 9 12 to 2% level is c omp arable in accuracy to that of the 
Cabibbo angle, which is about 1.4% at this moment |41[. Therefore, SADO can open the 
new era in which we can enjoy balanced knowledge of lepton and quark mixings. 

In the rest of this section, we focus on points which may have greater impact on wider 
areas of research, including symmetry tests, solar astrophysics and the earth science. Let us 
start by discussing implications to solar astrophysics, and then move on to the other topics. 
We assume, apart from Sec lVIIll D. the CPT invariance in our discussions. 

A. Observational solar astrophysics 

One of the purposes of the future solar neutrino experiments is to probe the deep interior 
of the Sun. We in fact live in an ideal location in the cosmos for this purpose, that is very 
close to the Sun, and the detailed information we can get for this one of the most standard 
main sequence stars should benefit wide area of astrophysics. In particular, observation of 
the full spectrum of the solar neutrinos which extends from 10 keV to 10 MeV region must 
shed light on our understanding of stellar evolution. It should provide us with the way of 
doing precise solar core diagnostics which is quite complementary to helioseismology. 

It is expected that pp neutrino flux will be calculated with error less than 2% because 
of the powerful luminosity constraint. Now what is the most serious obstacle for doing 
stringent test for such an accurate prediction? How accurately we can measure the flux of 
pp neutrinos which are about to leave the solar core? Assuming that the systematic error 
of the measurement can be controlled to less than 1-2% level, the largest uncertainty comes 
from the error in sin 2 9 12 determination. Therefore, we need an independent measurement 
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of 6*i2 to do precise observational solar physics. 6 Using the the most accurate value of 6*12 
which can be obtained by SADO, and using the observed solar neutrino flux at the Earth, 
we can determine accurately the infant solar neutrino flux at the solar core before they start 
to oscillate. 



B. Geo-neutrinos 

As described in Sec. IIII B[ the Earth is expected to be a very rich source of low energy 
[E v < 3.4 MeV) P e , whose detection can be of prime geophysical interest since they can 
provide otherwise inaccessible information on the abundance of radioactive isotopes such 
as U, Th and K inside the Earth and thereby help to unravel the internal structure and 
dynamics of the planet jlH Ht| • 

It has been shown that in a few years with a relatively small amount of exposure, Kam- 
LAND will have collected enough data in the energy interval 0.9 < -E pr0 mpt/MeV < 2.6 to 
clearly establish the presence of geo-neutrinos [3^. However, the precise determination of 
geo-neutrino flux requires considerably longer exposure (6 kt-yr to determine U flux within 
10% uncertainty [37 1). One of the reasons for the difficulty is that in this energy range, 
the neutrino flux from reactors is dominant at KamLAND and the geo-neutrino signal is 
predicted to be less than 1/5 of the reactor neutrino background. Therefore, precise mea- 
surement of 612 is necessary to reliably subtract the background contribution from reactors 
in order to determine precisely the flux of geo-neutrinos. 

The high-precision determination of #12 by SADO explored in this paper is essentially 
unaffected by the geo-neutrino background, as shown in Sec|V](See Fig.|SJ). Then, we can use 
this information to subtract the dominant contribution from reactors for the energy range 
relevant for geo-neutrinos, 0.9 < -Eprompt/MeV < 2.6, at KamLAND. We note that SADO 
itself would not be adequate to measure the geo-neutrino flux, because the geo-neutrino flux 
contribution in the energy interval 0.9 < .Eprompt/MeV < 2.6 is expected to be at most ~ 
10% at SADO. It is due to the fact that it is so close to the Kashiwazaki-Kariwa NPP, from 
which it receives a dominant contribution of reactor antineutrinos. 



C. Determination of masses and mixing parameters in future neutrino experi- 
ments 

The precise measurement of #12 should help identify the unknown quantities and to im- 
prove the accuracy of determination in future neutrino experiments. We mention here only 
two of them, neutrinoless double beta decay and CP phase measurement in LBL experi- 
ments. Neutrinoless double beta decay is probably the most promising way of identifying 
the absolute scale of neutrino mass in laboratory experiments |42|. As is well known the 
constraint imposed by the experiments on either the mass scale or the other observables 
such as Majorana phases involves 6i 2 . See e.g., [43] and the references cited therein. Its 
precise knowledge, therefore, should help identifying these quantities. 



6 Even in a strategy by which #12 can be determined simultaneously with the various components of solar 
neutrino flux by global fit |l 1| , precise information of #12 should improve the accuracies of solar neutrino 
flux determination. 
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Discovery of Leptonic CP violation in LBL experiments is one of the challenging goals 
in future neutrino experiments. The CP violating piece of the appearance oscillation prob- 
abilities P(Vfj, — > v e ) and Pip^ — > u e ) is proportional to Am|i and sin2#i2. Since detecting 
CP violating effects requires enormous precision in the experiments, uncertainties in these 
relevant parameters could easily obscure the discovery of CP violation. Moreover, once de- 
tection of CP violation is done, it will motivate and facilitate actual determination of the 
value of the CP angle 5. It does require very precise determination of Am^ and sin 2 6* 12 as 
well as Am^i and sin 2 #23, though the latter are not the subject in this paper. Therefore, 
accurate measurement of these quantities is the prerequisite for determination of 5. 



D. Test of CPT symmetry 

CPT symmetry is one of the most fundamental symmetries which are respected in rela- 
tivistic quantum field theory. It implies that the masses and the mixing angles of particles 
and their antiparticles are identical. Then, one can perform a stringent test of CPT symme- 
try in lepton sector by accurately measuring 9±2 and Arn^j with use of reactor antineutrinos 
and by comparing them with the same parameters measured by solar neutrino experiments. 

In the previous section we have shown that the sensitivity attainable in reactor measure- 
ment of 9 12 improves, in a reasonable setting, the one achievable by future solar neutrino 
experiments. Then, the significance of the CPT test is controlled by the accuracy reached by 
solar neutrino measurement. By comparing our predictions with those presented in Ref. [14J 
we can estimate the future sensitivity to CPT test in this sector. We note that the sensitivity 
on 9 12 in w]| essentially comes from solar neutrino data, not from KamLAND, as the latter 



is more efficient in decreasing the range of Am 2 2 but not of 9u ■ 

The current bound on the difference between sin 2 9 12 for neutrino and sin 2 9 12 for antineu- 
trinos is rather weak [H. Even if we assume that 9 12 is in the first octant, 7 



|sin 2 12 -sin 2 12 | < 0.3, (11) 

at 99.73% CL, which can be im pro ved by a factor 5 in the future if the accuracy expected 
from future solar neutrino data [14| is achieved and SADO is realized. See Table El 
For mass squared differences, the current bound j3] is 



I Am 2 x - Am 2 x I < 1.1 x 10" 4 eV 2 (12) 

at 99.73% CL, where Am 21 is the mass squared difference of antineutrinos. We observe 
that SADO will not be able to make significant improvement of this bound beyond what 
can be reached by combining future solar with KamLAND data since the bound will be 
essentially determined by the uncertainty of Am 2 2 coming from solar neutrino data which 
is significantly larger than that of Am| x determined by KamLAND or SADO. 

For comparison, let us mention that the only available constraint on # 2 3 is obtained by 
Gonzalez-Garcia et al. [45[ by analyzing atmospheric neutrino data. The bound obtained 
for the difference between neutrino and antineutrino mixing angles is: 

- 0.41 < sin 2 # 23 - sin 2 9 23 < 0.45, (13) 



7 If we allow 6*12 to be in the entire quadrant, < 6*12 < 7r/2, an another solution 8' 12 = tt/2 — 6*12 emerges 
Q|, and SADO's measurement will not be able to improve the CPT bound. 
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at 99.73% CL level. 

For Am§ 2 there are two results on CPT test, one by Super-Kamiokande group 
- 1.9 x KT 2 eV 2 < |Am 2 2 | - |Am| 2 | < 4.8 x eV 2 , 
and the other by Gonzalez- Garcia et al. |45[ 

- 10" 2 eV 2 < Am 2 2 - Am 2 32 < 3.4 x 10~ 3 eV 2 , 
both at 99.73% CL level 8 . 

E. Test of Non-standard neutrino interactions 

Precise measurement of Am^ and #i 2 by reactor experiments have an important impact 
in constraining the non-standard neutrino interactions (NSNI) with matter j42|. As dis- 
cussed recently in j48[, the currently allowed parameter region by solar and KamLAND 
data could be significantly modified were the NSNI present. It can be tested to a better 
accuracy by future data from KamLAND. The possible presence of such NSNI can be fur- 
ther constrained (or confirmed) by performing the precise measurement of Am 21 and 9 12 
with a SADO type reactor neutrino experiment. The point is that while solar neutrinos are 
severely affected by the matter effect reactor neutrinos are not. Therefore, such NSNI can 
be strongly constrained, assuming the CPT symmetry, if the mixing parameters inferred 
from solar and reactor neutrino observations coincide with each other. 



IX. CONCLUSION 

In this paper, we have investigated the potential power of dedicated reactor neutrino 
experiments for precision measurement of # 12 . By placing a detector in an appropriate 
baseline distance from a powerful nuclear reactor complex, and assuming 4% systematic 
error, a world-record sensitivity on sin 2 9\ 2 — 2% (~ 3%) at 68.27% CL is shown to be 
attainable by 60 GW t h-kt-yr (20 GWth-kt-yr) operation, superseding all the other proposed 
methods. Thus, it improves, after 20 GW t h-kt-yr operation the accuracy to be achieved 
jointly by KamLAND and the existing solar neutrino experiments for 3 years by more than 
a factor of 2. At 60 GW t h-kt-yr operation its enormous sensitivity of approximately 2% 
is about a factor of 2 better than that to be reached by additional 7 Be and pp neutrino 
observation with extremely small total errors of 5% and 1%, respectively (see Table HT|) . 

Toward the conclusion, we have carried out a careful estimation of the optimal baseline 
distance and obtained a rather wide range, between 50 and 70 km, from the reactor neutrino 



The bound by Super-Kamiokande was obtained by using the two-flavor analysis of atmospheric neutrino 
data, so that we have interpreted it to be the one placed on absolute values of Am| 2 . The bounds 
obtained in Ref. |45j are based on the analysis assuming the normal mass ordering for both neutrinos and 
antineutrinos. But, since both solar and KamLAND data imply that the splittings are hierarchical for 
both neutrinos and antineutrinos, and the 13 mixing angles for neutrinos and antineutrinos are favoured 
to be small, the bound is not expected to be quantitatively very different if neutrinos and/or antineutrinos 
masses have inverted ordering. 
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source for the current best fit value of Am^. The distance is determined by the require- 
ments that the first oscillation maximum occurs around a peak energy of the event number 
distribution in the absence of oscillation, and the geo-neutrino background is harmless. We 
have checked by taking a detector placed at 54 km from a reactor neutrino source that 
geo-neutrino background does indeed not produce any significant additional errors. 

To estimate the effect of background caused by other surrounding reactors we have exam- 
ined a concrete setting. We took the Kashiwazaki-Kariwa NPP in Japan, the most powerful 
reactor complex in the world as a reference, and assumed a detector (SADO) located in 
Mt. Komagatake which is 54 km away from the reactor complex. We have verified that 
the uncertainty due to the other 15 NPPs produce only 20% increase in the error of sin 2 #12 
determination (See Tab. [I}. We emphasize that since nuclear reactors are less populated in 
most of the rest of the world, examination of SADO with 15 NPPs will give us a conservative 
estimate of sensitivities for the similar reactor # 12 experiments on the globe. 

Of course, 8 12 is one of the fundamental parameters of particle physics and the signifi- 
cance of its precise determination by itself cannot be overemphasized. Nonetheless, we have 
discussed that there is a plethora of interesting physics implications available when such an 
accurate measurement of 6*12 is carried out. We point out its impact to solar astrophysics, 
geophysics, determination of the yet unknown mixing parameters, CPT symmetry test, and 
exploring non-standard interactions. 
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APPENDIX A: THERMAL POWERS AND THE DISTANCES TO THE DETEC- 
TORS 

In Table ITVl we present the thermal powers of the 16 NPPs in Japan and their distances 
to KamLAND and SADO at Mt. Komagatake which are used in our analysis. 

TABLE IV: The maximal thermal powers of the 16 NPPs in Japan and their distances to Kam- 
LAND and SADO at Mt. Komagatake are presented in units of GW t h and of km, respectively. 
Their fractional contributions of v e flux without oscillation to each detector, denoted as -Fkl.sadOj 
are also shown in %. The thermal powers of current operation are used for KamLAND analysis 
and the numbers and the distances are taken from For SADO sensitivity analysis, the future 
values for thermal powers are adopted for Hamaoka and Shika (indicated in parentheses), and the 
values are taken from the following web sites: 
http: / /www. chuden.co.jp/hamaoka/DETAIL/newgo-setsubi. html 
http : / / www . r ikuden .co.jp/ shika / outline2 / 

-^sado are calculated by using 1/25000 map, assuming that the Earth is a perfect sphere with 
radius of 6380 km. 



NPP 


Thermal Power (GWth) 


L KL (km) 


Fkl (%) 


^SADO (km) 


^SADO (%) 


Kashiwazaki 


24.3 


160 


32.0 


54 


77.0 


Ohi 


13.7 


179 


14.4 


354 


1.0 


Takahama 


10.2 


191 


9.4 


348 


0.8 


Hamaoka (future) 


10.6 (14.5) 


214 


7.8 


300 


1.5 


Tsuruga 


4.5 


138 


7.8 


299 


0.5 


Shika (future) 


1.6 (5.5) 


88 


8.2 


158 


2.0 


Mihama 


4.9 


146 


7.8 


307 


0.5 


Fukushima-1 


14.2 


349 


3.9 


144 


6.3 


Fukushima-2 


13.2 


345 


3.7 


139 


6.3 


Tokai-II 


3.3 


295 


1.3 


109 


2.6 


Shimane 


3.8 


401 


0.8 


492 


0.1 


Ikata 


6.0 


561 


0.6 


697 


0.1 


Genkai 


10.1 


755 


0.6 


828 


0.1 


Onagawa 


6.5 


431 


1.2 


236 


1.07 


Tomari 


3.3 


783 


0.2 


606 


0.08 


Sendai 


5.3 


830 


0.3 


953 


0.05 
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APPENDIX B: DETAILS OF ANALYSIS PROCEDURE 



Here we describe some details of our analysis procedure. We compute the expected 
number of v e events in the i-th energy bin, TV* 1160 = N* eac + Nf eo , where j\T reac and Nf eo are 
computed as follows 

iV7 eac (sin 2 # 12 , Am 2 21 ) = N p T 

/%pp i> 
dE v tk<Pk{E u )P{v e -> z/ e ; L k , E v )a{E u ) / dEe det R(E, E'){B1) 
k=i Ji 



and 



Nf °(sin 2 12 , <fo) = N p T (1 - X - sin 2 26 l2 ) 



x 



J dE v 1.83<k\fv(E v ) + fTt(E v )]a(E u ) JdEe det R(E,E'), (B2) 



where N p is the number of target protons in the detector fiducial volume, T is the exposure 
time, and <pk{E v ) is the neutrino flux spectrum from the fc-th NPP expected at its maximal 
thermal power operation. e& denotes the averaged operation efficiency of the fc-th NPP for 
a given exposure period and it is taken to be 100% here under the understanding that the 
unit we use GW t h-kt-yr refers the actual thermal power generated, not the maximal value. 
Pkipe —* v e ,Lk,E u ) is the familiar antineutrino survival probability in vacuum (given by 
Eq.(JU) with 6*i3 = and by setting /j = 8^) for the k-th NPP, and it explicitly depends on 
Am|i and sin 2 6*12. cr(E u ) is the v e absorption cross-section on proton, edet = 0.898 is the 
detector efficiency and R(E, E') is the energy resolution function, which is assumed to have 
Gaussian form, E (= E vxomp t) the observed prompt energy (total e + energy) and E' — E v - 
0.8 MeV the true one. 

We assume, following KamLAND 0, the energy resolution to be a(E)/E = 
6.2%/ a/ -E'(MeV) and consider each energy bin to be 0.425 MeV wide. The summation 
over k is meant to sum over the contributions from all the reactor NPPs, i.e., N NPP = 1 
or 16 depending upon whether only the Kashiwazaki-Kariwa or all 16 reactors are consid- 
ered. We have used in our calculations -Eprompt > 0.9 MeV and the future thermal powers of 
Hamaoka and Shika (see Table HV)) . Other relevant informations can be found in Ref. [igj ]. 

We have performed the geo-neutrino calculation as in Ref. 0u is the total flux of 

geo-neutrinos from U decays, we assume that the total flux of geo-neutrinos from Th decay 
is 83% of fa, f\j and /xh are the normalized energy distributions of U and Th geo-neutrinos. 
Here we have assumed the oscillation probability to be averaged. The procedure is not 
exactly correct but for our current purposes it gives a good approximation. 
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